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Abstract 

We obtain the equation of state of interacting alpha matter and the critical temper- 
ature of Bose-Einstein condensation of alpha particles within an effective scalar field 
theory. We start from a non-relativistic model of uniform alpha matter interacting 
with attractive two-body and repulsive three-body potentials and reformulate this 
model as a 0{2) symmetric scalar <p e field theory with negative quartic and positive 
sextic interactions. Upon restricting the Matsubara sums, near the temperature of 
Bose-Einstein condensation, to the zeroth order modes we further obtain an effec- 
tive classical theory in three spatial dimensions. The phase diagram of the alpha 
matter is obtained from simulations of this effective field theory on a lattice using 
local Monte-Carlo algorithms. 
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1 Introduction 



The a particle model of nuclei has its origin in the Gamow's theory of nuclear 
a emission, which assumes pre-formed a particles inside a nucleus. Light self- 
conjugate 4iV (N — 1,2, . . .) nuclei, called also a nuclei, have been associated 
with the a particle model, which assumes that a nuclei are composed of struc- 
tureless rigid a particles as their stable constituents; a itself is an extremely 
stable nucleus, and is emitted by naturally radioactive nuclei and the a nu- 
clei which are more stable than their neighbors [1]. This model in its entirety 
or some of its basic assumptions has been used over years in the a cluster 
structure investigations [2]. The a particle models are more successful in ex- 
plaining the excited states of a nuclei in terms of their rotational states than 
their ground state properties; a clustering appears as one of the many tran- 
sient types of clusterings that are continuously emerging and disappearing 
inside a nucleus [3]. 

The large N limit of a particle model describes a matter - an extended uniform 
system of structureless, indistinguishable a particles in a uniform background 
of neutralizing charge. Alphas maintain their integrity at densities where the 
inter-alpha distance is much larger than the average size of an alpha particle 
(< 2 fm) and at temperatures not too large compared to the binding energy 
per particle. Deconfined alpha matter occurs in astrophysical settings. At low 
densities/temperatures, characteristic for the nucleosynthesis site in supernova 
debris, the neutron to proton ratio and the ultimate outcome of the r-process 
nucleosynthesis depends on the recombination kinetics of nucleons into a par- 
ticles [4]. At higher densities (but still below the nuclear saturation density) 
a substantial fraction of hot and dense supernova matter resides in a parti- 
cles [5,6,7,8]. The equation of state of subnuclear matter [6,8] is at the heart 
of the simulations of the supernova collapse and explosions and has important 
implications for the formation of the supernova neutrino signal. 

The studies of a matter isolated from the astrophysical environment it is 
embedded in provide a first approximation for accessing the role of correla- 
tions, phase transitions and condensation phenomena within this component. 
There have been many detailed studies of the ground state of uniform a mat- 
ter [9,10,11,12,13,14,15,16,17], including the benchmark variational calcula- 
tions of Clark and co-workers [10,11,12,13] to which we compare our results. 
Alpha matter in its ground state may form a bcc lattice, which takes an advan- 
tage of closest possible packing of particles without destroying their identity 
[14,16,17]. Recent work focused on the possibility of the Bose-Einstein con- 
densation of a particles in a nuclei [18,19,20] and the competition between 
two- and four-body clustering in infinite matter [21]. 

The aim of this work is the study of uniform a matter from the lattice field 
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theory point of view. First, we show that the symmetric (non-condensate) 
phase of a matter interacting with attractive two-body and repulsive three- 
body forces, belongs to the universality class of 0(2) symmetric scalar 6 
field theories. Second, we simulate the scalar 6 field theory on a lattice using 
Monte-Carlo algorithms. The parameters of the continuum action are deduced 
from the non-relativistic interactions derived from the fits to the a -a scat- 
tering. For several parameter sets we then obtain the equation of state of 
a matter. Our approach is complementary to the other methods to study the 
a matter; an advantage of the present approach is the close connection to the 
field of critical phenomena, where many universal properties which are inde- 
pendent of details of interactions in the system (e.g. critical exponents) are 
known. 

Our work is in part motivated by the studies of the scalar 4 field theory which 
describe the shift in the critical temperature of the Bose-Einstein condensa- 
tion for repulsive atomic systems [29,30,31,32,33,34]. These studies are suited 
for treating repulsive, dilute Bose gases in the vicinity of the critical temper- 
ature of Bose Einstein condensation. Alpha matter is attractive at large and 
repulsive at small distances. To describe these features, we adopt the scalar 
6 theory with negative quartic and positive sextic couplings. While we con- 
centrate here on alpha matter we anticipate that our model could be applied 
to dilute Bose gases interacting with long-range attractive and short range 
repulsive interactions. Section II derives the effective action of the classical 
three-dimensional theory starting from the Hamiltonian of a -matter inter- 
acting with two- and three-body forces. In Section III we describe the Monte- 
Carlo simulations of this theory on a lattice. The results for the equation of 
state of a -matter and the critical temperature of Bose-Einstein condensation 
are presented in Section IV. Section V is a short summary of this work. 



2 Effective action 



A uniform, non-relativistic a matter is described by the Hamiltonian 



H 



2m 



V^ t (x)V^(x) - /i^ f (x)^(x) + ^ f (x)^(x)f/(: 



where m is the a particle mass, /i is the chemical potential, ^(x) and ^(x) 
are the creation and annihilation operators, and 



3 



U(x) = J dxV 2 (x', x)^t(xXx') 

+ Jdx!j dx'V 3 (x, x', x")^ t (x")^(x')^ t (x")V(x'), (2) 

where V% and V3 are the two- and three-body interaction potentials (which are 
not identical to the free-space interactions in general, see below). To reformu- 
late the theory above as an effective field theory, we replace the finite range 
potentials by zero range constant couplings (which are constrained to repro- 
duce the low-energy scattering data) and adopt the lattice regularization. For 
contact form of interaction Eq. (2) reduces to 

*7(x) = <fe^t( x )^(x) +^ t (x)V(x)^ t (x)V(x), (3) 
where gi = Anh 2 a^} jm and the two-body scattering length is defined as 



For systems featuring short-range repulsive interactions (hard core) the two- 
body potential in Eq. (4) should be replaced by the scattering T-matrix which 
sums up the particle-particle ladder series. We follow the convention that 
assigns postive sign to the scattering length for attractive interactions. In the 
case, where the three-body interactions are represented by a sum of pairwise 
(repulsive) two-body interactions: V 3 (x, x', x") = (x — x') + V^ 1 (x — x") + 
V^^x' — x") [23], the three-body coupling constant is #3 = 47ih 2 a^) /3m, 
where a^) is defined via a spatial integral (4) over the two-body potential 

vf (x-*0- 

To map the theory defined by Eq. (1) on an effective scalar field theory, we 
adopt the finite temperature Matsubara formalism, described in refs. [29,35]. 
Consider the fields ip and ip^ as periodic functions of the imaginary time 
r G [—/3,j3], where (3 = 1/T is the inverse temperature. These fields can 
be decomposed into discrete Fourier series 



i!>(x,u> v )= £ e^(x,r) 

v=— OO 

00 

= Vo(x)+ £ e^(x,r), (5) 

u=—oo, uj^O 

where the Fourier frequencies u v are the bosonic Matsubara modes u v = 
lixivT {y assumes integer values). The Matsubara Green's function for our 
system can be written as 
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G M (u u , x) = [wj v - (2m)~ 1 V 2 + (6) 

where the chemical potential includes contribution from the self-energy. To 
clarify this point, consider the spectrum of excitations of the system which 
is determined above the critical temperature T c by the poles of the retarded 
propagator. This function is obtained from the Matsubara Greens function via 
the replacement iu v — > uj + irj (where uj is real) 



G R (u, x) = [u - (2m)" 1 V 2 - Re£ fl (x, u)+n f + 177] -1 , (7) 

where /i/ is the chemical potential of non-interacting gas and E ii (x, a;) is 
the retarded self-energy. The contribution of the two-body interaction to the 
retarded self-energy can be written as 



Ef 2) (x,u) = jrfxWr^t^x'.u + w'jlmG^x'w'), (8) 

where Fu) is the four-point vertex, which to the lowest order reduces to the 
potential V2, but in general may include two particle resummation series. The 
three-body contribution to the self-energy, which results from the three-body 
repulsive interaction, can be written in terms of six-point vertex rVg) in similar 
fashion. For contact interactions T( 4 ) = g 2 and IVg) = g 3 and the net self-energy 
is a constant S . Thus the retarded propagator can be written as 



G r (uj,x) = [uj - (2m) _1 V 2 + fi + irj]- 1 , \i = y. f - E . (9) 

For temperatures T > T c a macroscopic number of particles starts to occupy 
the zero-momentum state. Since the chemical potential must stay below the 
lowest possible energy state to keep the occupation probabilities positive, it 
follows from Eq. (9) that fi < 0; because a macroscopic number of particles 
occupies the zero-momentum ground state as T — > T c we need fi 0. Since 
the above argument is equivalent to the statement that the fugacity 



z = e? 11 -> 1, (10) 

the characteristic scales of spatial variations of the Green's function with non- 
zero Matsubara frequencies are / = (2muj v ) ~ 1 ^ 2 , which are of the order of 
the thermal wave-length A = (27r/mT) 1 / 2 . The contribution of the non-zero 
Matsubara modes to the sum in Eq. (5) can be neglected if one is interested 
in the scales L ^> I, which are characterized by the zero-frequency modes. 

Upon dropping the second term in Eq. (5) and introducing two new real scalar 
fields (pi and (f>2 
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where 77 = Jm/h 2 f3, the continuum action takes the form 



(11) 



/^UE[^0(x)] 2 



r / / \2 u r / / \2i 2 w r 1 1 \ 
+ 2 0(x) " 4! i 0(x) J + 6! ^ 0(x) 




(12) 



where <p 2 = 4>\ + 0|, r = -2/3/J7/ 2 , u = 4!/3g 2 ?7 4 , if = 6!/3^ 3 ?7 6 . The action (12) 
describes a classical 0(2) symmetric scalar <fi 6 field theory in three spatial 
dimensions. Note that the positive sextic interaction guarantees that the en- 
ergy is bound from below, which otherwise would be unbound because of the 
negative sign of the quartic term that describes the attractive two-body inter- 
actions. The characteristic length scale of the theory is set by the parameter 
u, which has dimension of inverse length; the dimensionless parameter of the 
lattice theory is ucll, where is the lattice spacing. The theory described by 
the action (12) differs from the ordinary 4 theory by the order of the phase 
transition it predicts. Indeed, the 6 theory admits a first order phase tran- 
sition since at low temperatures the free energy develops a second minimum 
with the field </> assuming two equilibrium values, one at the origin (<p) = 



and a second along the circle of constant radiusy (0 2 ) + (4%). While the first 
order (liquid-gas) phase transition in a matter is of interest in its own right, 
below we concentrate on a the possibility of the Bose-Einstein condensation 
of a particles. This phenomenon is characterized by a pole in the distribu- 
tion of Bose particles, which then occupy the zero momentum ground state 
of the system - a state that is characterized by the condition \x — > at fixed 
temperature and density. 

The thermodynamic functions of the model are obtained from the partition 
function 



For example, the expectation value of the particle number density n a = 
(i[)*ip) = (/3V)~ 1 dlriZ/dfi, where V is the volume. 

The coupling constants of the continuum theory can be related to the values 
of the two-body scattering length a sc (derived from a -a potentials fitted 
to the scattering data) only in the dilute limit. The details of short range 
physics become unimportant when the condition A ^> a sc is fulfilled, where 
A 2 = 27r / ?7 2 is the thermal wavelength of a particles. However, the concept 





(13) 
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of scattering length can be used for dense systems if it is defined in terms 
of the full scattering T-matrix, instead of the Born amplitude (see e.g. [26]). 
Since the two-body a -a potentials are characterized by hard cores (strong, 
short-distance repulsion) we shall treat the scattering length as a parameter 
with the understanding that it can be derived from an appropriate T-matrix 
(or phase-shifts), where the effects of the hard-core are removed by a full 
resummation of the ladder series. 



3 Numerical simulations 

For the purpose of numerical simulations we discretize the continuum theory 
on a lattice by replacing the integrations over a summation over the lattice 
sites. The discretized version of the continuum action (12) is 



where the z-summation is over the sites of a 3d cubic lattice, the v summation is 
over unit vectors in three spatial direction (nearest neighbor summation). The 
hopping parameter k, the two- and three-body coupling constants A and £ are 
related to the parameters of the continuum action through the relations a\r = 
(1 — 2X)/k — 6, A = aLK 2 u/6 and ( = wk 3 /90. The lattice and continuum fields 
are related by 0l( x ) = (2K/a^) 1//2 0(x). The components of the spatial vector 
x v assume integer numbers of the lattice spacing a^. x M = 0, a^, . . . (L u — l)a^. 
In our numerical simulations we adopt a simple cubic lattice in three spatial 
dimensions with periodic boundary conditions imposed on the field variable: 
^(x + olL) = <Pl{x) [for a box of length L there are L 3 (real) variables within 
the volume (Lai,) 3 ]. 

The lattice regularization of a field theory provides an effective cut-off A ~ aj} 
and breaks its translational and rotational invariance. To restore the imparl- 
ance and recover the continuum theory the limit ua^ — > should be taken in 
the numerical data. In the context of 4 theory in four (and higher) spacetime 
dimensions it has been proven that the continuum limit of the theory is trivial, 
i.e. the renormalized couplings vanish in the continuum limit (this would not 
be inconsistent for an asymptotically free theory). In three dimensions the 4 
and 4> 6 theories are superrenormalizable. We leave the study of the renormal- 
ization group equations governing the present 6 theory for the future. Here 
we take the point of view of the effective field theory, which exploits the fact 




2^0 L (x)0 L (x + ai>)0 L (x) 2 + A 1 + 0l(x) 




(14) 
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that the cut-off can be made so large (spacing so small) that the cut-off effect 
in the Green's functions of the theory are completely negligible [24]. In other 
words, we require the correlation length £ of the theory to satisfy the condition 
£ ^> ci£. The correlation length is deduced by computing the correlation func- 
tions on the lattice (<f>(x)4>(y)) for a fixed lattice spacing a^. In the simulations 
described below it was found that fixing the physical lattice spacing cll = 1 
fm guarantees that the condition £ ^> cll is fulfilled; this would correspond 
to introducing a hard core of a radius ~ 1 fm in the interactions, which is 
consistent with the size of the cores of a -a potentials. 



Another limit that must be taken in the lattice Monte-Carlo simulations is 
the limit uL — > oo to minimize the finite lattice size effects; in other words, 
the thermodynamic parameters of the system in a finite box must have a limit 
when the size of the box goes to infinity, and this limit should not depend on 
the periodic boundary conditions imposed. In practice we successively increase 
the size of the lattice and compare the observables (e.g. the expectation values 
of the fields and the correlations functions) until a convergence is reached. The 
finite size effects become insignificant when the condition L ^> £ is fulfilled. 



Within the effective field theory approach sketched above the "bare" param- 
eters r, u, and w can be used to derive the parameters of the lattice action 
for a specific choice of lattice spacing a L . We note here that the parameters 
of the continuum action themselves are renormalized due to the higher order 
Matsubara modes. These modifications, which are expected to be small in the 
vicinity of the critical temperature of phase transition, are neglected. 



A universal feature of the potential models of the nuclear a -a interaction is 
a repulsion for small separations, which is weakening as the relative orbital 
angular momentum / increases from to 4, followed by an attractive tail that 
is largely independent of I. The range of the inner repulsion is approximately 
2 fm, the range of the outer attraction is about 5 fm. For the hard-core plus 
a square well potential of Van der Spuy and Pienaar (VP) [25], with the core 
radius a = 1.7 fm, the well depth Vo = 7.2 MeV and well width 4 fm, one finds 
the S wave scattering length a sc = —7.33 fm which translates into a coupling 
constant g 2 = —954.3 MeV fm 3 . Below we shall vary the magnitude of a sc 
to account for the renormalization by the short range effects and differences 
among the available a -a interactions. For the three body coupling constant 
we employ the potential of Ref. [23] to obtain g% = 5944.8 MeV fm 3 . Table 1 
lists the parameters of the continuum and the lattice theories which were used 
in the numerical simulations for = 1 fm. 
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0.01 0.02 0.03 0.04 

n [fm 3 ] 

Fig. 1. Dependence of the energy (per particle) of the a matter on the density 
for several values of the scattering length at constant fugacity log z = —0.1. The 
squares and the triangles show the variational results of Johnson and Clark [13] for 
the potentials of ref. [25] and [27], respectively. 

4 Results 

The field configurations on the lattice were evolved using a combination of 
the heatbath and local Metropolis algorithms with the number of equilibra- 
tion sweeps 10 5 — 10 6 . Lattices with sizes from 8 3 to 64 3 were used. The results 
shown below were obtained for the lattice size 32 3 (the changes in the observ- 
ables, as the lattice size is further increased, are insignificant). Once the field 
values on the lattice were determined, these were transformed into their coun- 
terparts in the continuum theory to obtain the statistical average value of the 
Hamiltonian (H), i.e. the grand canonical (thermodynamical) potential Q as 
a function of density. The energy of the a matter per particle Q a /N a is shown 
in Fig. 1 for the three values of the scattering length a^' = —7.33, -10, and -12 
fm and the same three-body interaction. For comparison we show also the zero 
temperature variational results of Ref. [13]. As expected, the thermodynami- 
cal potential increases with the scattering length, since the two-body potential 
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Fig. 2. The density dependence of the critical temperature of a condensation for 
scattering lengths a sc = 12 fm (triangles) and o sc = 7 fm (squares). The critical 
temperature for isoscalar quasideuteron pairing [28] is shown by dashed dotted line. 

scales linearly with a sc . The close agreement of our result for a sc = —7.33 fm 
derived from the VP-potential with the results of Johnson and Clark obtained 
for the same interaction suggests that at low temperatures neither the temper- 
ature dependence of the thermodynamic potential Q a /N a , nor the contribution 
from / > partial wave are significant (taking the limit T — > in the present 
theory is complicated by the appearance of the Bose-Einstein condensate). 
Our result for the larger scattering length a sc = — 12 fm fit quite well the 
results of Johnson and Clark based on the Ali-Bodmer (AB) potential [27]. 

The critical temperature T c of the Bose-Einstein condensation of a particles 
can be obtained from the simulations. In practice we obtain the function n a (T) 
rather than directly T{n a ) at constant fugacity z — > 1. Since we work at small 
log z = —0.1, we compute the density of the system at constant small chem- 
ical potential and temperature which is identified with T c . In this manner 
we extract an approximate value of the critical temperature as a function of 
a matter density. Fig. 2 shows the density dependence of the critical tempera- 
ture of Bose-Einstein condensation of a matter extracted from simulations for 
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Table 1 

The parameters of the continuum and lattice theories for several values of the tem- 
perature and two-body scattering lengths at constant log z = —0.1. The parameters 



of the lattice theory are listed for 


cll = 1 fm. 






T 


r [fm 2 ] 


n/10 3 [fm- 1 ] 


w/10 5 


K 


A 


c 


= -7.33 fm 


9.0 


0.174 


0.959 


3.113 


0.047 


-0.355 


0.361 


8.0 


0.154 


0.853 


2.460 


0.049 


-0.348 


0.331 


7.0 


0.135 


0.746 


1.883 


0.052 


-0.340 


0.299 


6.0 


0.116 


0.640 


1.384 


0.056 


-0.330 


0.265 


5.0 


0.096 


0.533 


0.961 


0.060 


-0.318 


0.228 


ag) = -10.0 fm 


9.0 


0.174 


1.309 


3.113 


0.041 


-0.372 


0.244 


8.0 


0.154 


1.164 


2.460 


0.043 


-0.366 


0.224 


7.0 


0.135 


1.018 


1.883 


0.046 


-0.359 


0.204 


6.0 


0.116 


0.873 


1.384 


0.049 


-0.350 


0.181 


5.0 


0.096 


0.727 


0.961 


0.053 


-0.339 


0.158 


5.00 


0.10 


0.73 


0.96 


0.053 


-0.339 


0.158 


a|c } = -12.0 fm 


9.0 


0.174 


1.571 


3.113 


0.038 


-0.382 


0.193 


8.0 


0.154 


1.397 


2.460 


0.040 


-0.376 


0.178 


7.0 


0.135 


1.222 


1.883 


0.043 


-0.369 


0.162 


6.0 


0.116 


1.048 


1.384 


0.045 


-0.361 


0.145 


5.0 


0.096 


0.873 


0.961 


0.049 


-0.350 


0.126 



two values of the scattering length. These temperatures are above the critical 
temperature of the pairing in the 3 Si- 3 Di pairing in a non-interacting nucleon 
matter [28]. While at low densities the a condensation is the dominant effect, 
it should be kept in mind that the a particles lose their identity at higher 
densities when their wave-functions strongly overlap and the Pauli-principle 
starts to work for the constituent nucleons. The critical density for the ex- 
tinction of a particles is sensitive to the magnitude of the spatial extension 
of their wave-functions (n* oc R~ 3 ) and is of the order of n* ~ 0.03 fm -3 
for an effective radius R a ~ 2 fm. At densities above the so-called Mott den- 
sity a particles are destroyed due to the Pauli principle acting among their 
constituent nucleons [21]. 
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5 Summary 



The finite temperature theory of a matter close to the critical temperature of 
the Bose-Einstein condensation has been reformulated as a scalar Euclidean 
(p 6 theory with a negative quartic and positive sextic interactions. The theory 
is restricted to the zeroth order Matsubara modes and is described by a real 
two-component field which is 0(2) symmetric in the order parameter space. 
The Monte-Carlo simulations of the theory on the lattice established numer- 
ical relations between the thermodynamic parameters (we work at constant 
fugacity z « 1 and extract the density, the energy density and the critical 
temperature). The simulations are carried out within an effective field theory 
with a cut-off, which satisfies the condition <C £, where £ is the correlation 
length. The obtained equation of state is in a good agreement with the previ- 
ous variational calculations of the binding energy of infinite a matter at zero 
temperature [10,11,12,13] in the intermediate density regime. The equation of 
state of low density alpha matter is presented here. 

The critical temperature of Bose-Einstein condensation of a particles is found 
larger than the temperature for BCS condensation in the isoscalar 3 Si- 3 Di 
(quasi-deutron) channel. The dominance of a condensation is a low density 
phenomenon, since at densities above the so-called Mott density a particles 
are destroyed. The mechanism is the Pauli principle acting among their con- 
stituent nucleons [21]. A treatment of the Mott effect requires the knowledge 
of the internal structure of the alpha particles and is beyond the scope of 
present formalism. Thus at asymptotically large densities (of the order of nu- 
clear saturation density) one is left with a quasi-deutron BCS condensate [28]. 

This work concentrated on a system consisting exclusively of a particles. How- 
ever, the Hamiltonian (1), which contains two-body attractive interaction, may 
have a spectrum of bound states. It is known that 12 C (the three alpha state) 
is bound, therefore our system can evolve into a gas of 12 C particles, and even 
further to the point where the matter consists of 56 Fe - the most stable nu- 
cleus. As pointed out in the Introduction, calculations of the matter properties 
under supernova conditions show a substantial fraction of alphas present in 
the matter, which motivates the present study. Nevertheless, it will be inter- 
esting to study the spectrum of bound states of the Hamiltonian (1) and, for 
example, the peculiarities related to the spectrum of the three-body bound 
state (Efimov effect) [36]. 

We are grateful to Prof. J. W. Clark for communicating to us refs. [12,13] and 
to Dr. Kurt Langfeld for discussions. This work was in part supported by a 
grant from the SFB 382 of the DFG. 
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